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Abstract 

In the present study we revise the possible polaron contribution to the charge and energy 
transfer over long distances in biomolecules like DNA. The harmonic and the simple inharmonic 
(U(x) = x 2 /2 — /3x 3 /3) lattices are considered. The systems of PDEs are derived in the continuum 
approximation. The PDEs have the one-soliton solution for polarons on the harmonic lattice. It 
describes a moving polaron, the polaron velocity lies in the region from zero to the sound velocity 
and depends on the polaron amplitude. The PDEs describing polarons on the inharmonic lattice 
also have the one-soliton solution only in the case of special relation between parameters (parameter 
of inharmonicity (3 and parameter of electron-phonon interaction a). Polaron dynamics is numeri- 
cally investigated in the wide range of parameters, where the analytical solutions are not available. 
Supersonic polarons are observed on inharmonic lattice with high inharmonicity. There is the range 
of parameters a and (3 where exists a family of unusual stable moving polarons with the envelope 
consisting of several peaks (polarobreather solution). The results are in qualitative agreement with 
recent experiments on the charge transport in DNA. 



Keywords: polaron, charge transport, DNA 
PACS numbers: 71.38.-k; 87.15.-v 



I. INTRODUCTION 



This paper continues our investigations on the polarons in ID lattices [1]. Its primary 
goal is an explanation of the high efficient charge transfer in biological macromolecules. 
Charge transfer (CT) is of utmost importance in physics, chemistry and biology (2 1 3 1 . 
or instance, it is essential in p rocesses at nano- and micro-scales, for solar cells technology 
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molecular electronics 



15 



16] and other novel technologies. The CT is very important, 



e.g. for oxidative DNA degradation 17]. Noteworthy is that DNA has been labelled as 
electric conductor [18], semiconductorjl9 1 . insulator J20I and even superconductor 21 |. 
Such wide range of possible electric properties is explained by the variety of different DNA 
structures used, diversity of experimental methods and inner and outer surrounding DNA 
conditions. It has been shown that the rate of CT exponentially decays with the distance L 

22^24|. 



as kcT oc exp(— 7L), where 7 > 0.1 — 0.2 A 



-1 



Significant progress in enhancing the conducting properties of DNA was achieved after the 
DNA synthesis reached its present status of automation. Apparently adenine (A) is the most 
^referable building block because of its resistance to charge trapping j^j] and its CT efficiency 
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27). Three sets of duplexes were synthesized giving 4-, 6- and 14-base A tracts with a 



covalently attached rhodium complex [Rh(phi)2(bpy')] 3+ serving as the photooxidant [28]. It 
was found that there was no significant change in degree of decomposition (7 = 0.0013 A -1 ), 
i.e. the efficiency of charge transport very slowly depends on the distance, what contrasts 
with the larger values found in previous experiments. The possibility of CT over 34 nm 
(100-mer DNA) was recently announced by J. Barton and collaborators [29] and this result 
surpasses most previous achievements with molecular wires. 

Analogous results were obtained for the synthetic q-helical peptides where a very shallow 



30] . One interesting feature in 



distance dependence of the charge transfer was also found 
these studies was revealed: the charge transport is a coherent single-step process, i.e. the 
charge carrier can be transferred ballistically 31 



Traditionally the hopping and tunnelling are considered as a most abundant mechanisms 



in the theory of charge transfer in DNA 
candidates for charge carriers 
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-[36]. Polarons were also extensively studied as 



44j (see also review |45|). At present, the polaron theory 



of CT in biomacromolecules is very popular. Polarons are considered in different models 
and approximations 
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Effects of nonlinearity can give contribution to the charge transport due to formation 



of the bounded soliton-electron state 
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55] or polarobreathers 56N64]. As a working 



definition breathers are localized solutions with at least two degrees of freedom in a one- 
dimensional lattice (i.e., they have an "internal" degree of freedom, a vibrational motion) 



57j. 



while stable localized solutions with one degree of freedom are solitons 

Recently the concept of nonlinearity was additionally advanced and solectrons, - bounded 
state of solitons and charge carriers, were found in numerical experiments (see 65] and 
references therein). 

Notably, the charge transfer in DNA appearers to be independent of distance but is 
critically sensitive to different kinds of defects, - static and dynamical disorder 66|, l67| . 



Biomacromolecules are very complex systems and many kinds of fluctuations can disturb 
the path-ways of charge migration. A DNA macromolecule is supposed to have a uniform 
or a random distribution of the helix angles. And an influence of angles fluctuations on the 
charge transfer efficiency was analyzed 
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71]. Dynamical fluctuations (temperature) also 
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can affect the hopping probability thus giving rise to charge localization 

In the present paper we consider the charge transport in harmonic and inharmonic lattices 
and derive analytical expressions for moving polarons in the continuum approximation. 
Polaron dynamics is investigated numerically in the range of parameters where the analytical 
solutions are not available. Special attention is given to possible explanation of recent 
experiments on the charge transport in DNA. 



II. MOVING POLARON ON THE HARMONIC LATTICE 

We start from consideration of ID harmonic lattice with free boundaries consisting of N 
particles with the hamiltonian 

N , N-1 

3=1 3=1 

—t 

where \I/ = ipi, -0 2 > ■ ■ ■ > V'iv is the discrete wave function, m is the particle mass, k is the lattice 
rigidity, Xj is deviation of the jth particle from the equilibrium. "Particle" represents the 
DNA base, and their interaction is due to the tt—tt overlapping of neighboring heteroaromatic 
bases. 
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The components of the electron-phonon interaction operator are: 

H tj = ( \j:i< j + Sjj-itj-i + Sjj+itj, (2) 

where 5ij is the Kronecker symbol; tj - hopping integral between jth and (j + l)th particles; 
Cj is the energy of interaction of the charge carrier with the lattice site (on-site energy). 
This strategy is known as the tight-binding approach to DNA which takes into account the 
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on-site energies and suitably parametrized hopping onto the neighboring sites 

The main object of the present investigation is synthetic DNA comprised of identical base 
pairs. So the diagonal disorder is absent and ej = const. This constant is the electronic 
energy origin and without the loss of generality can be put to zero. The hopping integral is 
often written in the form suggested by Su, Schrieffer and Heeger 
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tj = -[v - a(x j+1 - xj)] , (3) 

where vq is the hopping integral at the equilibrium and parameter a accounts the electron- 
phonon interaction. The electron charge transfer is considered, though all results are valid 
for the hall transfer. 

It is convenient to make hamiltonian ([T]) dimensionless using mass m, energy Vq and the 
rigidity coefficient k as the units. Then the unit of time [t] = y^m/k, unit of length [L] = 



y/vo/k. If the DNA parameters are chosen 38, 39], namely vq = 0.3 eV, k = 0.85 eV/A 2 , 
m = 130 a.m.u., then the time unit is (t) ~ 0.13 ps, and the length unit is [L] ~ 0.59 A. 
Parameter a is also dimensionless: a/\/v k — > a. The dimensionless value of this parameter 
for DNA a ~ 1.2. Parameter a is the single parameter specifying the lattice. 
The dimensionless hamiltonian (CQ) reads: 

^ N ^ N-l N-l 

3=1 3=1 3 = 1 

where qj = (xj+i — Xj) is the relative displacement of neighboring particles. This hamilto- 
nial together with the Schrodinger equation iM> = H e ^> generates the system of evolution 
equations: 

%j =(lj - Qj-i) ~ oc[(ip*ipj+i + c.c.) - (ipj^j + c.c.)] 

• i (5) 

^3 =~ K 1 - acQj-im-i + (1 - uqjWj+i] , 
h 

~ h Fk ~ 

where h is the dimensionless Planck's constant: — \ > h f» 0.017. 

v o V m 



The evolution equations in variables qj and ipj are 

Qi ={Qj+i ~ 2 1j + Qj-i) ~ + cc.) - 20*^+i + c.c.) + (ipj^ipj + c.c.)] 

• i (6) 

h 

Eqs. (jBJ) are more useful for analytical consideration, while Eqs. are more suitable for 
numerical calculations. 

Eqs. (JSJ) can be reduced to the system of continuous equations in the long-wave approx- 
imation. Let us expand the variables gj and ipj into Taylor series: 

i I I 2 1 // i 

Qj±i = Qj ± e?j + e^±- 

1 , (7) 
^i±i = i>j ± <^ + -e 2 ip" ± . . . 

where e is small parameter. After substitution of Eqs. ([7]) into Eqs.flU]) one can get the 
system of nonlinear partial differential equations (PDEs): 

qtt =Qxx + 2a(il)ip*) xx 

(8) 

ifH =~[2(l-ag)V + Vw], 
h 

where terms of the order e 3 and higher are omitted. This system is not integrable, but it 
has special one-soliton solution: 

.4 




cosh 2 [d(x - u p t)] 
B exp[i(kx + ut)] ^ ' 



cosh[<i(x — v p t)] 

where 1/d is the polaron width and v p is the polaron velocity. A, B are amplitude of relative 
displacements and amplitude of the wave function, correspondingly, (kx + ut) is the phase of 
;he wave function. If v p = then Q) coincides with the solution for the unmovable polaron 
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Solution is the one-parametric solution with the following relation between parame- 
ters: 

/ — 7 / fo?\ 1 fo? „ [d . . 

Wt " 1 -2\Ia ; B = Vr (10) 



and 



~2 £ < 1 ; w = — ~ — »!■ (n) 




a) b) 

FIG. 1: (a) Dependence of the unmovable polaron amplitude A vs. parameter a. Solid line 
is analytical prediction, dots are results of numerical simulation, (b) Dependence of the polaron 
velocity v vs. the amplitude A for two values of parameter a. Analytical results (|10p are shown by 
solid (a = 0.4) and dashed (a = 0.7) lines. Numerical results are shown by circles (a = 0.4) and 
squares (a = 0.7). 



and the amplitude of relative displacements A can be chosen as a free parameter specifying 
the polaron. Other parameters (d, v p ) are defined through amplitude A. The norm of the 

oo 

wave function J \ip(x)\ 2 dx = 1 was preserved while deriving (TiT)]) - ffl~Tj) . 

— oo 

The polaron velocity v p lies in the range from zero (unmovable polaron) to the sound 
velocity (the dimensionless sound velocity t> sn d = 1). The amplitude of the moving polaron 
A > A min and A min = a 3 is the amplitude of unmovable polaron. 

Now we check the analytical results in numerical simulation of Eqs. ([5]). First of all we 
find the range of parameter a where the continuous approximation is valid. Fig. [1^ shows 
the dependence of the unmovable polaron amplitude vs. parameter a of the electron-phonon 
interaction. Analytical and numerical results coincide for a < 0.4. 

Fig. [T|d shows the dependence of polaron velocity vs. its amplitude A for two values 
of parameter a. As expected, coincidence is very good for a = 0.4 and there is some 
disagreement for larger values of a. But despite some discrepancy for a = 0.7, the numerical 
data form a distinct dependence. It means that there can exists the polaron solution of ([6]) 
which, nevertheless, is not governed by the continuum approximation. 
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a) b) c) 

FIG. 2: Snapshots of the polaron evolution on the harmonic lattice at different time moments: 
t = (a), t = 250 (b) and t = 500 (c). Polaron is centered on the site jo = 25 at t = 0. Positive 
values along the ordinate axis are the modulus of wave function (solid line), and negative values 
are relative displacements (dashed line). Lattice parameters are a = 0.4, N = 500. Initial polaron 
amplitude A = 0.4. 



We use parameter a = 0.4 for the analysis of the polaron stability. Figs. [2] demonstrates 
high polaron stability, where snapshots of the polaron evolution are shown at three time 
moments. The initial conditions are chosen according to (|9|)- (TTTj) for A = 0.4. The polaron 
is very stable, it travels « 400 lattice sites without any noticeable change. The calculated 
polaron velocity (v p ~ 0.78) is in good agreement with analytical prediction (i> aria i = 0.775). 

The polaron initially generated according to the "correct" initial conditions is very stable. 
The problem is whether the polaron can be formed and is it stable in the case of arbitrary 
initial conditions? To answer this question the initial excitation is chosen according to 
(15])- (111)) for A = 0.4, but the velocity is v p = 0.2 instead of correct value t> ana i = 0.775. 
The evolution of this initial state is shown in Fig. |3j Stable polaron with parameters A rs 
0.21, v p ~ 0.68 forms soon (polarons in Figs. [3b and [3b are practically identical). The 
polaron amplitude decreased relative to initial value. But relations f fTUl) and ([Til are fulfilled 
for parameters of this polaron with high accuracy, i.e. the analytical polaron velocity v axia \ ~ 
0.67 for the polaron amplitude A = 0.21. The initial excitation transforms into polaron 
accompanied by the noise irradiation. Noteworthy, that 100% of the electron wave function 
is localized at polaron; no fraction of the wave function is irradiated with the noise. 

Wide range of initial conditions was tested and in all cases stable polaron is formed with 
the relations between parameters obeying (110p - (llip . 
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a) b) c) 

FIG. 3: Polaron evolution with initial conditions: a = 0.4, A = 0.4 and polaron velocity v p = 0.2 
which does not satisfy correct relations (|lUp . Snapshots are shown at t = 0, t = 250 and t = 500. 
Lines and axes are the same as in Fig. [2j N = 500. 



In the next section we transform Eqs. (|S]) into the integrable system of PDEs. 



(12) 



79] where a special variable 



III. REDUCTION TO THE INTEGRABLE SYSTEM OF PARTIAL 
DIFFERENTIAL EQUATIONS. 

Let us rewrite system (jSJ) for convenience: 

qtt =q xx + 2a(ipip*) xx 

tpt =~[2(1 - aq) ip + ip xx ) . 
h 

We follow the reductive perturbation method (RPM) 
transformation 

q = eui + e u 2 + ... 

y =s 1/2 (x - Ct) 

r=e 3 /\ 

is made and e is small parameter. The main objective of the RPM usage is the reduction of 
the degree of derivative by time q a in ( TT2l) to the first order. The new (laboratory) coordinate 
system moves with velocity c relative to the old coordinate system (second line in (fl~3l) ). 

Substituting (fl~3l) into (fl2"l) and equating terms of the order e 2 one gets that c = 1. It 
means that the laboratory coordinate system moves with the sound velocity c = f sn( j = 1 
relative to the old coordinate system. 



(13) 
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Equating terms of the order e 3 and taking into account that the solution tends to zero at 
±00, we get 

' 2u T = 2a(ipip*) y 

(14) 

ih(i/j T — IPy) = —2(1 — atl)tjj — tjjyy. 

Additional linear variables transformation (ip(y,T) = a^(z,T) exp[i(kz — wT)], u(y,r) = 
bU(z,T), t = cT, y = dz) is necessary to reduce (IT4"|) to the "standard" form: 

U T =2(^^*) z 

(15) 

where coefficients labelled with 'tilde' are some linear combinations of "old" coefficients. 

System (TT5]) is exactly solvable and its one-, two-soliton etc. solutions are known (80]. As 
a particular case, the one-soliton solution has the same functional form as solution fl9]). But 
the relation between parameters A, B, d, v p is somewhat different. Numerical simulation 
confirms high stability of polarons described by the solution of ( Tl5l) . 

The numerical experiments of polaron interaction shows that the polarons preserve their 
shapes after collision. This fact additionally proves the fact that integrable system ([15]) is a 
good approximation for discreet system ([6]). 

In the next section we consider the moving polaron on the inharmonic lattice. 



IV. MOVING POLARON ON THE INHARMONIC LATTICE 

In this section we consider the inharmonic lattice with the dimensionless FPU potential: 

U(q) = \q 2 ~ (16) 

where /3 is the parameter of non-linearity and q is the relative displacement of particles. Many 
realistic potentials (Morse, Lennard- Jones, etc) are reduced to f fl~6]) when the deviations from 
equilibrium are not large (the expansion of potentials into the Tayler series up to the third 
order). Moreover, the a-FPU potential (Tl6]) allows to make necessary analytical calculus. 

We employ the same procedure as in the case of harmonic potential. We expand the 
relative displacements and the wave function into Taylor series according to ([7]). The result 
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is the system of PDEs: 

Qtt = (qxx + j^QxxxxJ ~ (3 (q 2 ) xx + 2a ; 



(17) 

i/; t =~[2(1 - aq) - i) xx }. 
h 



Unfortunately, the relation between coefficients in ffTTl) is such, that this system can not be 
reduced to the exactly solvable PDEs [80] which has the canonical form: 

3u tt ={u xxx - Quu x ) x + 8\4>\l x ; 

(18) 

i 4>t =4>xx - u<p. 

However, if (3 = 2 a (coefficients of the second spatial derivatives in the first equation in 
( |T7|) are equal), then system (fT7|) has the special one-soliton solution. This solution has the 
same form as for the harmonic lattice: 

.4 



q(x,t) 
ip(x,t) 



cosh 2 [d(x — vt)} 
B exp[i(kx + out)} 



(19) 



cosh[<i(a; — vt)] 
with the relations between parameters: 



d = y/aA = y/WA, 



(5A\ 1 ' 2 1 1 [tf (5 A 




B = 
and 

hv v , 2 + d 2 -k 2 

k = — £ <1; oo = = >1. (21) 

2 h 

There is an addition term in the expression for the polaron velocity (I20p as compared to the 
polaron velocity on the harmonic lattice (fTUj) . This positive term is due to the nonlinearity, 
and the actual velocity is determined by the balance of electron-phonone interaction and 
nonlinearity. In the absence of the electron-phonon interaction (a = 0) the velocity is 



v ~ 1 + (3 A/6, and is the soliton velocity on the a-FPU lattice 81 1. 

Now we verify the polaron stability in numerical simulations on the inharmonic lattice in 
a way analogous to the harmonic lattice. Fig. Ufa shows the dependence of the unmovable 
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FIG. 4: (a) The dependence of the polaron amplitude vs. parameter a for the unmovable polaron. 
Solid line is analytical solution (|19p with parameters (|20p -([2i p . dots are numerical results, (b) 
Dependence of the polaron velocity vs. amplitude A for two values of parameter a: a = 0.4 
(solid line is analytical dependence, circles are numerical results), a = 0.7 (dashed line is analytical 
dependence, squares are numerical results). 

polaron amplitude vs. the parameter a of the electron-phonon interaction. One can see, 
that there is also a good agreement between analytical and numerical results for a < 0.4. 

Figure Hb shows the dependence of the polaron velocity vs. amplitude A of relative dis- 
placements for two value of parameter of electron-phonon interaction a. A good agreement 
between analytical and numerical results is also observed for a = 0.4 but there is some 
gap for a = 0.7. Similar to the harmonic lattice, the numerical data form a well-defined 
dependence for a = 0.7. It means that there can exist polarons which are not governed by 
the continuum approximation. 

Numerical simulations show high polaron stability. Figs. |5] illustrates snapshots of the 
polaron evolution at three time moments, when the initial conditions are generated according 
to (fl~9l) -( l21~|) . The initial polaron amplitude is A = 0.3. The polaron preserves initial 
parameters in the course of the evolution. The calculated polaron velocity is v p w 0.78, 
whereas the analytical velocity t> ana i — 0.76 for A = 0.3. 

The polaron evolution with improper initial conditions is shown in Figs. |6j The initial 
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FIG. 5: Snapshots of the polaron evolution on inharmonic lattice (a = 0.4, f3 = 0.2) at different 
time moments: t = (a), t = 250 (b) and t = 500 (c). The polaron is centered on the site jo = 25 
at t = 0, the initial polaron amplitude A = 0.3. Positive values along the ordinate axes are modulus 
of the wave function \ipj\ (solid line), negative values are relative displacements qj (dashed line). 
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a) b) c) 

FIG. 6: Polaron evolution with improper initial conditions. Snapshots are shown at t = (a), 
t = 250 (b) and t = 500 (c) . Notation and axes are the same as in Fig. [5] 

polaron is formed according to Eqs.( ll9p -( l2~Tj) . but the polaron velocity is chosen to be v p = 0.2 
instead of analytically predicted t> ana i = 0.76. The simulation show that the polaron self- 
organizes (compare polarons in Figs. |6h> and [6fc). Parameters of self-organized polaron 
satisfy f l20]) -( l2~Tj) with high accuracy. The "numerical" polaron has velocity and amplitude 
v p ~ 0.65 and A w 0.16 whereas t> ana i = 0.63 if A = 0.16 for the "analytical" polaron. 

The numerical simulations show that polaron is stable on the inharmonic lattice. How- 
ever, the results demonstrated above are valid only if a — 2 /3 and a < 0.4. In this case 
the lattice inharmonicity is small ((3 = 0.2) and does not influences essentially on the po- 
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laron evolution. The case of large a and f3 are of greater interest. This case is considered 
numerically in the next Section. 



V. NUMERICAL SIMULATION ON THE NONLINEAR LATTICE 
WITH ARBITRARY PARAMETERS 

Equations 

xj =(?j - <?j-i) + /3(gJ_i - gj) - a[(ip*tij j+1 + c.c.) - + C - C -)L 

• < (22) 
^ == [(1 - aqj-iji/jj-i + (1 - a?j)^-+i] ; g,- = - Xj 
h 

are integrated numerically with different initial conditions. We use initial conditions in the 
form f fT9|) . but A, d, v are arbitrary parameters and k = 0. The question is whether polarons 
are really exist in the case when the continuous approximation is not valid. 

Fig. [7J illustrates how the strange four-peaked polaron-like excitation self-organizes on 
the inharmonic lattice with parameters a = 0.4, = 1.0. The polaronic nature of this 
excitation is supported by the 100% bearing of the wave function by the localized lattice 
compression. Initial excitation at t = is chosen with parameters d = 0.1, A = 0.6, v p = 0.7, 
k = 0. Lattice excitations behind the polaron do not have the electronic components and 



are solitons. The linear relation between soliton velocities and amplitudes 81] is clearly 
visible in Fig. [7] and additionally supports the solitonic nature of these excitations. This 
numerical simulation is an example of how an arbitrary initial conditions transform to the 
polaron with extra perturbations emitted in the forms of polaron. Noteworthy that there 
are no fluctuating lattice vibrations. 

To verify the stability of the obtained four-peaked polaron, we use it parameters 
(xj, Vj, ipj] j = 980 -T- 1020), located between two horizontal arrows in Fig. [7J as the ini- 
tial condition for new numerical simulation. The result of evolution is shown in Fig. [HJ The 
polaron is very stable. After travelling through w 400 lattice sites its parameters are not 
changed and it moves with constant supersonic velocity v p ~ 1.14. 

The most natural and expected shape of the polaron envelope is a smooth bell-shaped 
form and an existence of four-peaked polaron is quite unexpected. But there exist polarons 
with other number of peaks. 

Fig. [9] shows how the two-peaked polaron self-organizes. The initial excitation is chosen 
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FIG. 7: The result of the evolution of the initial excitation at t = 800. The formed polaron is 
indicated by two arrows for clarity. There are six solitons behind the polaron. The front of the 
sound propagation is shown by the vertical arrow. The initial excitation is centered at jq = 50. 
Modulus of the wave function - solid line; relative lattice displacements - dashed line. TV = 1200. 
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FIG. 8: The evolution of the four-peaked polaron formed from the initial conditions of previous 
numerical experiment. Inserts show polarons with higher resolution. N = 500. 



more narrow than in the previous numerical experiment. The initial parameters are d = 0.5, 
A = 0.6, v p = 0.7, k — 0. The evolution of this initial state is more complex: there 
are intermediate one-, two- and three-peaked excitations. And finally, after rather long 
time, the two-peaked polaron becomes the stable state. Its velocity is somewhat less then 
the sound velocity. Numerical simulations show that other initial conditions can produce 
three-peaked polarons. 

And finally we consider the inharmonic lattice with parameters a and /3 corresponding 
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FIG. 9: The snapshots of the polaron at different time moments: t = 321 (a), t = 345 (b),i = 380 
(c), t = 400 (d). 



to the parameters in DNA. Usually the Morse potential is used to account the interaction of 
neighboring bases U(q) = D[l — exp(— bq)] 2 . An expansion of this potential into the Taylor 
series gives the a-FPU potential. In the dimensionless form the potential is U (q) ~ ^ 3 
with f3 1.1 and a ^ 1.2. 

This set of parameters does not allow an analytical consideration. Numerical simulations 
are carried out for wide range of initial conditions. As an examples the evolution with 
different initial parameters are analyzed. Results are shown in Fig. [TUJ Polarons self- 
organize from all used initial conditions. But no multi-peaked polarons are found at these 
parameters choice. 
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FIG. 10: The result of evolution of the excitation with different initial parameters at t = 350: 
A = 0.6, d = 0.4, v p = 1.0 (a); A = 0.6, d = 0.1, v p = 1.2 (b); A = 0.5, d = 0.2, v p = 0.5 (c). 
Lattice parameters: N = 500, a = 1.2, (3 = 1.1. 



The obtained results are rather unusual. First, the polaron can be supersonic with velocity 
exceeding the fastest soliton velocity. Second, the polaron envelope can be multi-peaked and 
consists of two, three and four peaks (multi-peaked polarons were found earlier in 84| in 
a slightly different problem formulation). It means that the minimum energy states can 
have a breather (multi-peaked) or a solitonic (single-peaked) character and depend on the 
parameter values and initial conditions. Third, the wave function is fully concentrated inside 
the self-organized polaronic potential well, and emitted solitons do not carry away a minute 
amount of a wave function. 

Polarons on the homogeneous lattice were considered above. But DNA and polypeptides 
can have static and dynamical defects of different nature. An interaction of polarons with 
lattice defects is considered in the next section. 



VI. CONCLUSIONS 

In conclusion we briefly summarize the main results. 

The moving polarons on the harmonic and inharmonic lattices are considered. An ana- 
lytical solutions for the moving polarons are derived in the continuous approximation. This 
approximation is valid for the large-radius polarons and when the electron-phonon inter- 
action a is weak and the nonlinearity parameter is small. In few limiting cases these 
equations are reduced to the exactly solvable models with special one-soliton solutions. The 
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polaron on the harmonic lattice is subsonic. The polaron on the inharmonic lattice can 
travel with the supersonic velocity. An excellent agreement between analytical results and 
numerical simulation was obtained. 

In the intermediate range of parameters a and /3, where the continuous approximation is 
not valid, the family of stable solutions with unexpected properties was found in numerical 
simulations. These solutions are supersonic with velocities exceeding the velocity of the 
fastest solitons. The envelop of these polarons for the relative displacements consists of 
few peaks. These polarons are similar to polarobreathers as they have internal vibrational 
structure. For the larger parameter values, typical for DNA, all formed polarons are subsonic 
with the bell-shaped envelope. 

In all cases the wave function is concentrated with the 100% probability inside the polaron 
potential well. 

The results can be useful in explaining the recent experiments on the highly efficient 
charge transfer in synthetic oligonucleotides, where the CT occurs as the single-step coherent 
process. 
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